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Abstract. Excitation functions A]si{piab, ©cm.) of the analyzing power in elastic proton-proton scattering have been 
measured in an internal target experiment at the Cooler Synchrotron COSY with an unpolarized proton beam and a po- 
larized atomic hydrogen target. Data were taken continuously during the acceleration and deceleration for proton kinetic 
energies Tiat (momenta pjab) between 0.45 and 2.5 GeV (1.0 and 3.3 GeV/c) and scattering angles 30° < 0c.m. < 90°. 
The results provide excitation functions and angular distributions of high precision and internal consistency. The data 
can be used as calibration standard between 0.45 and 2.5 GeV. They have significant impact on phase shift solutions, in 
particular on the spin triplet phase shifts between 1.0 and 1.8 GeV. 



PACS. 25.40.Cm Elastic proton scattering - 13.75.Cs Nucleon-nucleon interactions 
nomena in reactions - 21.30.-x Nuclear forces 



24.70.-I-S Polarization phe- 



1 Introduction 

' The internal target experiment EDDA (3 01 at the Cooler 
' Synchrotron COSY 0] is designed to provide high precision 
measurements of proton-proton elastic scattering excitation func- 
tions ranging from 0.45 to 2.5 GeV of laboratory kinetic energy 
' Tiab- In phase 1 of the experiment, spin-averaged differential 
cross sections were measured continuously during acceleration 
with an internal polypropylene (CH2) fiber target, taking par- 
ticular care to monitor luminosity as a function of beam mo- 
mentum 1 1, 3]. In phase 2, excitation functions of the analyz- 
ing power Aff are measured using a polarized atomic hydrogen 
beam as target. In phase 3, excitation functions of the polariza- 
tion correlation parameters Aatjv, Ass and Asl are measured 
using a polarized beam and a polarized target. The present pa- 
per gives a detailed account of our measurements of the ana- 
lyzing power excitation functions Ajv(p(a6, 0c.m.) using an in- 
ternal polarized atomic beam target during beam acceleration 
and decelaration, see also |2]. 

Following the Argonne notation we denote the analyzing 
power by An |5] where the spin directions are defined with 
respect to the laboratory frame of reference {S, N, L) with N 
the normal to the scattering plane, L the longitudinal (beam) 
direction and S the sideways direction (S ^ N x L). In the 
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Saclay notation |a] the laboratory frame of reference is denoted 
(s, n, k) with (s, n, k)^{S, N, L). The analyzing power mea- 
sured with beam polarization normal to the scattering plane is 
denoted {N, 0; 0, 0)=AooiiO- If the target is polarized the ana- 
lyzing power is denoted (0, N; 0, O)=yloooji- In PP elastic scat- 
tering beam and target particles are identical and Aoono=^ooon 
=An. Often, the symbol A or P is simply used for A^. Since 
AN{0c.m.) is antisymmetric with respect to ©cm. =90°, the 
An data are presented only in the forward cm. angle range 
0°-90°. 

Elastic PP scattering experiments ff\ are fundamental to the 
understanding of the NN interaction. For kinetic energies be- 
low 1 GeV a precise database of differential cross sections and 
polarization observables has been accumulated. These data are 
well represented by phase shift solutions |8, 9, 10, 11, 12, 13, 
14]. Modern phenomenological and meson-theoretical poten- 
tial models LI 5^ .16. .17.. .ISt. .191 provide excellent descriptions 
of the data up to the pion threshold. Extending the meson- 
theoretical models to higher energies requires the inclusion of 
inelastic channel contributions. Using relativistic transition po- 
tentials and restricting to the NN, NA and AA channels yields 
reasonable descriptions of the data up to about 1 GeV I Tsll. 
Those models can be improved by including other nucleon res- 
onances than the A. But at certain higher energies the meson 
exchange model has to break down when the hadron substruc- 
ture reveals itself in a crucial way i20ll2Tll . Besides meson ex- 
change models recent theoretical work is based on chiral per- 
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turbation theory (xPT) IH pi pi MM 113 HE IH HI, 

Skyrmion or Soliton models I3ll l32l l33 l and quark model 
descriptions |35, 36, 37, 38,S|43Ell|4S, 43, 44, 45] of the 
NN interaction. A recent review of the theoretical progress can 
be found in l-^. 

Elastic pp scattering at GeV energies is ideally suited to 
study the short range part of the NN interaction. At 2.5 GeV 
kinetic energy a four momentum transfer of up to 1.5 GeV/c is 
reached corresponding to spatial resolutions of about 0.13 fm. 
The precise knowledge of the analyzing power as a function of 
energy provides a focus on heavy meson exchanges, especially 
on the role of the cj-meson with respect to the spin-orbit force. 
Apart from the true meson-exchange genuinely new processes 
might occur at small distances involving the dynamics of the 
quark-gluon constituents. Another issue related to the quark- 
gluon dynamics is the question of existence or nonexistence 
of dibaryons. Various QCD inspired models predict dibaryonic 
resonances with cm. resonance energies Eji ranging between 
2.1 and 2.9 GeV. Not any resonance has been observed so far. 

A complete listing of previous analyzing power measure- 
ments in the kinetic energy range 0.45-2.5 GeV can be found 
in the SAID database L14il . Recent measurements at discrete ki- 
netic energies Tiah>l GeV are from the SATURNE II faciUty 
f47, 48, 49, 50, 51, 52]. At lower kinetic energies Tiab<1.0 GeV 
high precision measurements of the analyzing power at discrete 
energies have been performed at LAMPF and re- 

cendy at lUCF |5l|5l|5l|5l. 

A first attempt to measure excitation functions of the an- 
alyzing power using an internal target during beam accelera- 
tion has been performed at KEK |60, 61]. However, data were 
taken only at one fixed lab-angle, Oiat — 68°, from 0.5 to 
2.0 GeV using a polarized beam and an unpolarized target. In 
this experiment two narrow structures have been observed near 
Tiab^632 MeV (p,ah=1259 MeV/c), i.e. in the neighbourhood 
of a depolarizing imperfection resonance with 7G=3 of the 
KEK ring. These structures were not confirmed in an external 
target experiment at SATURNE II 

The motivation of the EDDA experiment was to measure 
excitation functions at intermediate energies in a large angular 
range with a high relative accuracy. As a first result EDDA pro- 
vided (j| excitation functions da/df2{piab, Oc.m.) of unpolar- 
ized pp scattering. Addition of these internally consistent data 
to the SAID database 1 12] allowed to extend the global PSA 
from 1.6 GeV to 2.5 GeV jlj. 

The technique applied by EDDA to measure excitation func- 
tions during the acceleration of a proton beam is perfectly ap- 
propriate to provide a new polarization standard in the form of 
precise excitation functions of the analyzing power Aat. In ad- 
dition, the measurement of analyzing powers is a first necessary 
step in the EDDA program towards measuring spin correlation 
parameters Amn, Ass and Asl as a function of energy. With 
a polarized hydrogen target, a high and stable polarization is 
available and internally consistent analyzing power data can be 
taken over a wide energy range. In sect. II we describe the ex- 
perimental setup and the modifications of the EDDA detector 
to meet the increased demands for vertex reconstruction. The 
data acquisition and processing is presented in sect. III. The 
results are given in sect. IV and discussed in sect. V. 



2 Experimental Setup 

2.1 Overview 

The EDDA experiment is designed to provide high precision 
measurements of the proton-proton elastic scattering excitation 
functions over a wide energy range. Using an internal target, 
data taking proceeds during the synchrotron acceleration ramp 
of COSY, so that a complete excitation function is measured 
during each acceleration cycle. Statistical accuracy is obtained 
by averaging over many thousand cycles (multi-pass technique). 
This technique requires a very stable and reproducible opera- 
tion of COSY. The internal recirculating COSY beam provides 
beam intensities high enough for use of a polarized atomic 
beam target. Typical values are 3 • 10^" unpolarized protons 
in the ring, recirculation frequencies of 1.20 - 1.57 MHz and 
target densities of 2 • 10^^ hydrogen atoms/cm^ yielding lumi- 
nosities of 7.2 - 9.4 • lO^"^ cm^^g-i 

The analyzing power measurements were performed using 
an unpolarized proton beam and a polarized hydrogen target. 
This method is in contrast to the usual method of using a po- 
larized proton beam and an unpolarized hydrogen target. It has 
the advantage to avoid all uncertainties and systematic errors 
due to depolarization resonances in the acceleration ramp. The 
direction of the target polarization was changed from cycle to 
cycle between ±x and ±y thus allowing a proper spin flip cor- 
rection of false asymmetries ]63]. The absolute value of the 
target polarization is constant during the time period of an ac- 
celeration cycle (15 s). Small drifts of the absolute polarization 
over periods of up to two weeks are taken into account. There- 
fore, excitation functions can be measured with a high relative 
accuracy. 

2.2 Detector 

The EDDA detector shown in Figure [l] in a schematic fash- 
ion consists of two cylindrical detector shells. The solid an- 
gle coverage is 30° to 150° in Oc.m. for elastic proton-proton 
scattering and about 85 % of An. In phase 1 of the experiment 
only the outer detector shell was used to take unpolarized dif- 
ferential cross section data with a CH2 target (and C target for 
background subtraction) ]1, 3J. The outer detector shell i64ll 
consists of 32 scintillator bars (B) which are mounted paral- 
lel to the beam axis. They are surrounded by scintillator semi- 
rings (R) and semi-rings made of scintillating fibers (FR). The 
scintillator cross sections were designed so that each particle 
traversing the outer layers produces a position dependent sig- 
nal in two neighbouring bars and rings. The resulting polar and 
azimuthal angular resolutions are about 1° and 1.9° FWHM. 

For the measurement of spin observables, however, a polar- 
ized atomic beam target shall be used. The interaction region of 
such a target is far from being pointlike. The polarized hydro- 
gen beam has a finite diameter of about 12 mm (FWHM). It is 
superimposed on a background of unpolarized hydrogen atoms 
which are due to the residual gas in the beam pipe vacuum. 
Therefore, a second inner detector shell is required to provide 
an appropriate vertex reconstruction. 

The inner detector shell (HELIX) is a cylindrical hodoscope 
consisting of four layers of 2.5 mm diameter plastic scintil- 
lating fibers which are helically wound in opposing directions 
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Fig. 1. Scheme of the EDDA detector. The outer hodoscope consists 
of scintillator bars B, scintillator semi-rings R and semi-rings made of 
scintillating fibers FR. The inner hodoscope HELIX consists of four 
layers of scintillating fibers helically wound in opposing directions. 



Fig. 2. Scheme of the polarized atomic hydrogen beam target. 



so that coincidence of hits in the lefthanded and righthanded 
helices gives the point at which the ejectile traversed the ho- 
doscope. The 640 scintillating fibers are connected to 16-channel 
multianode photomultipliers and read out individually using 
the LeCroy proportional chamber operation system PCOS III. 
A detailed description of the helical fiber detector can be found 
in ll65ll . Its purpose is (i) vertex reconstruction of elastic proton- 
proton scattering events in conjunction with the outer shell, 
(ii) background suppression of scattering events from the back- 
ground of unpolarized hydrogen atoms surrounding the polar- 
ized hydrogen beam and (iii) improved angular resolution. Com- 
bined with the spatial resolution of the outer detector shell, the 
helix fiber detector provides for vertex reconstruction with a 
FWHM resolution of 1.3 mm in the x- and y-direction and 
0.9 mm in the z-direction. Using a fit of the vertex and scat- 
tering angles with constraints imposed by pp elastic scattering 
kinematics the resulting polar and azimuthal angular resolu- 
tions are about 0.3° and 1.3° FWHM. 

2.3 Polarized Atomic l-lydrogen Beam Target 

The polarized atomic hydrogen beam target l66ll is shown in 
Fig. 12] Hydrogen atoms with nuclear polarization are prepared 
in an atomic -beam source with dissociator, cooled nozzle, per- 
manent sixpole magnets and RF-transition unit. The design of 
the atomic beam target had to meet constraints imposed by the 
EDDA experiment. First, the space close to the interaction re- 
gion is limited: target components must be outside the angular 
acceptance of the EDDA detector This dictates a rather large 
distance of about 30 cm between the second sixpole magnet 
and the interaction region. Second, in view of the closed or- 
bit distortions in COSY only weak guide fields are allowed 
such that only one pure hyperfine state can be used. Third, the 
COSY beam width at injection and the change of the horizon- 
tal beam position during acceleration makes the use of a storage 
cell (typical apertures 10 mm x 10 mm) unfavorable. Reducing 
the beam width with beam scrapers, and thus the injected beam 
current would at least partly offset the benefit of higher target 
densities and lead to increased background. 

In the dissociator, hydrogen is dissociated in an inductively 
coupled 350 W RF-discharge and passes through an aluminum 



nozzle cooled to about 30 K and a skimmer. About 5 • 10 
hydrogen atoms per second are produced. For the differential 
pumping a turbomolecular pump of 600 liter/s is used in the 
first vacuum chamber between nozzle and skimmer and four 
turbomolecular pumps of 1000 liter/s in the subsequent vac- 
uum chambers. In addition, a cryopump with a pumping speed 
of about 2300 liter/s is installed in the second sixpole vacuum 
chamber. In the beam dump turbomolecular pumps of 360 and 
600 liter/s and a cryopump of 1000 liter/s are installed. The 
comparable low temperature of the nozzle leads to a decreased 
velocity (most probable velocity 1.3 km/s) of the atomic beam 
and thus an increased target thickness. The atomic beam tar- 
get is usually operated at 0.5 mbar liter/s hydrogen flow. A 
small amount of oxygen flow is mixed into the hydrogen flow 
yielding a thin layer of H20-molecules in the region of the 
cooled nozzle. Thus, the recombination of hydrogen atoms by 
the cooled nozzle is minimized resulting in an increased target 
density. 

The atomic beam source selects hydrogen atoms in a pure 
hyperfine state (mj=+l/2, m/=+l/2), in order to achieve a 
high polarization in a weak magnetic holding field. Here, nij 
and TO/ are the magnetic quantum numbers of the electron and 
proton spins, respectively. Hydrogen atoms in the TOj=+l/2 
state are focused by the sixpole magnets while those in the 
TOj=— 1/2 state are defocused. The first sixpole magnet re- 
moves the two hyperfine states with m j =— 1 /2 and the Abragam- 
Winter RF-transition unit induces an intermediate field transi- 
tion (IF-transition) to a depopulated hyperfine state, (TOj=+l/2, 
TO/=— 1/2)^ (mj=— 1/2, m/=+l/2), which is removed by 
the subsequent sixpole magnet. 

In Fig. |3l typical intensity and polarization distributions are 
shown as a function of the longitudinal position Zy. The plot 
shows the elastic proton-proton scattering rates to the left and 
right of the detector (a) and the derived polarization (b). The 
polarized atomic hydrogen beam stands out at z = on an 
unpolarized background which is due to the residual hydrogen 
gas in the beam pipe vacuum. Polarizations of about 90 % are 
deduced when the unpolarized background is subtracted. The 
atomic beam width (FWHM) is about 12 mm at the intersection 
with the COSY beam. A guide field in the interaction region of 
1 .0 mT pointing up or down was used to align the target spins in 
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Fig. 3. (a) Elastic proton-proton scattering rates to the left and right of 
the detector and (b) the derived polarization profile as a function of the 
longitudinal vertex position 2„ measured with vertical target polariza- 
tion at Tiab=793 MeV. The polarized atomic hydrogen beam stands 
out at Zv=0 on an unpolarized background of the residual hydrogen 
gas in the beam pipe vacuum. The dashed lines indicate the Zv vertex 
cut [-15,-1-20] mm. The decreasing number of events at z„< — 15 mm 
is due to the decreasing detector acceptance. 



the vertical direction. The decrease in rate below z<— 15 mm is 
due to the acceptance of the EDDA hardware trigger. Including 
the unpolarized background, the effective polarization is about 
73% for a z-vertex cut [-15 mm, +20 mm] and the effective tar- 
get thickness is 1.8-10^^ atoms/cm^. The background of unpo- 
larized hydrogen depends on the quality of the COSY vacuum. 
Therefore, in addition to the standard ion getter pumps, tita- 
nium sublimation pumps are switched on. Then, the resulting 
vacuum in the region of the EDDA-target is < 1 • 10^^ mbar 
without and about 1 • lO^'* mbar with hydrogen beam target. 

The atomic beam profile in y-direction has been measured 
at a fixed flattop energy (796 MeV) by sweeping the COSY 
beam by steerer magnets across the target and measuring the 
vertex distribution of elastic proton-proton scattering events. 
The target density distribution was deduced by fitting a Gaus- 
sian distribution for the polarized atomic beam plus a con- 
stant unpolarized hydrogen background. The finite width of the 
COSY beam and the detector resolution function were taken 
into account by an appropriate folding of the distributions. The 
resulting atomic beam width (FWHIVI) in y-direction of about 
12 mm compares well with the corresponding width in z-direc- 
tion. Similarly, the polarization profile in T/-direction which was 
deduced from the measured left-right asymmetries compares 
well with the corresponding profile in z-direction. 

Due to the background of unpolarized hydrogen atoms the 
effective polarization depends on the overlap of the COSY beam 
with the polarized atomic beam. Fortunately, variations of the 
vertical position and width of the COSY beam during the ac- 
celeration are very small, see Figs. [S] and |6l In addition, the 
COSY beam is much smaller than the polarized atomic beam. 
Therefore, the effective target polarization is practically con- 
stant during the acceleration and deceleration. 

In the beam dump the polarization of the atomic beam is 
continuously monitored using an additional RF-transition unit 
and a permanent sixpole magnet as Breit-Rabi polarimeter 




ferrite yoke 



Fig. 4. Magnet for the guide field of the target polarization. Magnetic 
guide fields of 1.0 mT in ±x- and ±j/-direction are produced near the 
beam axis by linear superposition of A- and B-fields. 



2.4 Magnetic Guide Field 

The magnet for the guide field of the target polarization is 
shown schematically in Fig.0] It consists of a ferrite yoke and 
ferrite pole shoes located external to the vacuum chamber The 
magnetic fields in x- and y-direction are produced by superpos- 
ing two fields of equal strength along the azimuthal directions 
45° and 135°. The earth magnetic field is shielded to a high 
degree by the ferrite yoke. The residual magnetic field is com- 
pensated by a small correction field using the magnet for the 
guide field. The currents for the magnetic guide fields and the 
field corrections are controlled by a personal computer The di- 
rection of the target polarization was changed from cycle to cy- 
cle by changing the direction of the guide field from +x to —x 
and +y to —y. Using a flux gate probe the resulting field distri- 
butions were measured as a function of x, y and z. These mea- 
surements proved a high magnetic field homogeneity across the 
fiducial vertex volume. 

The strength of the magnetic field was chosen to be 1.0 mT. 
This value is sufficiently large compared to the field distortion 
produced by the circulating beam particles at the interaction 
point. On the other hand, the resulting distortion of the closed 
orbit is sufficiently small. It can be calculated using the COSY 
lattice parameters at the target point. For beam momenta be- 
tween 1.0 and 3.3 GeV/c the resulting angle kick varies be- 
tween 60 and 18 /irad yielding a horizontal closed orbit shift 
Ax between 5 1 and 16 fim for guide fields along the ?/-direction. 
Similar values are obtained for the vertical closed orbit dis- 
tortions. The closed orbit distortions deduced from the vertex 
reconstruction of the data are in good agreement with these es- 
timates. 



2.5 COSY Beam 

The ramping speed of COSY was changed from the usual value 
1.1 (GeV/c)/s to 0.55 (GeV)/s and data were taken during ac- 
celeration as well as deceleration of the COSY beam. The time 
period of one cycle was about 15 s. With an average of 2.8-10^° 
protons in the ring, luminosities of about 8-lO^^cm^^s^^ were 
achieved and accumulated to an integrated luminosity of ^10^^ 
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cm^^. The beam parameters were continuously measured dur- 
ing the acceleration ramp. The beam momentum is derived 
from the RF of the cavity and the circumference of the closed ^ 
orbit with an uncertainty of 0.25 to 2.0 MeV/c for the lowest E 
and highest momentum, respectively. " 10000 
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Fig. 5. Vertex distributions in the {x, ?/)-plane for the running period 
November 98. 



The x-y density distributions of the COSY beam can be de- 
duced from the measured vertex distributions. Thus, beam po- 
sition and width can be reconstructed as a function of the beam 
momentum. Typical vertex distributions are shown in Fig.|5]in 
the form of two-dimensional scatter plots. Since the polarized 
atomic beam target is directed in x-direction and its width in 
j/-direction is rather large these distributions resemble approx- 
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Fig. 6. Beam position (thick lines) and beam width (FWHM, thin 
lines) in x- and y-direction as a function of the accelerator time for 
three running periods. The horizontal lines mark the start of data tak- 
ing in the acceleration ramp, the beginning of the flattop momentum 
region and the beginning of the deceleration ramp. 



imately the COSY beam density. Interestingly, the horizontal 
beam position and width change slowly during the accelera- 
tion and deceleration whereas the vertical beam position and 
width stay nearly constant. This behaviour can be seen more 
clearly by plotting the mean values and x, y and the FWHM- 
widths Ax—2.35(jj^ 
cycle time t, see Fig.|6l 

The deviations from the nominal beam position x = and 
?; = indicate (i) misalignments of the detector and (ii) closed 
orbit distortions of the COSY beam which are caused by small 
deviations from the nominal magnetic dipole fields in the ac- 
celerator ring. In all three running periods the horizontal devi- 
ations X vary systematically as a function of beam momentum 
between at most -5 mm and +2 mm. Vertically, a small but con- 
stant offset of y « +1.5 mm was observed, see Fig.|6l 

Variations of the horizontal and vertical beam widths 
and (jy are expected since the adiabatic damping causes a l/piab 
dependence of the beam emittances and Cy. In addition the 
optics of the accelerator ring, i.e. the amplitude functions /3x 
and (3y may depend on piab- The equations for cTj. and (jy read 



Ay=2.35(7y as a function of the COSY 



(1) 



If the amplitude functions ™d f3y are constant a 1/ y/piab 
dependence of the beam widths is expected. However, in the 
chosen mode of operation (3.j., was decreasing and (3y increasing 
with piab- As a consequence, the horizontal beam width was 
rather strongly decreasing with piab whereas the vertical beam 
width ay was nearly constant. 
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Fig. 7. Typical COSY machine cycle: beam momentum p=piab, beam 
current I and trigger rate as a function of the COSY cycle time. The 
dashed vertical lines denote the flattop region. 



The approximate constancy of the vertical beam profile is 
of great importance for the effective target polarization which 
depends on the overlap between beam and target. The vertical 
width of the target beam (12 mm, FWHM) is about a factor two 
larger than the vertical width of the COSY beam and the vari- 
ations of the vertical beam position and beam width during the 
acceleration and deceleration are negligibly small with respect 
to the vertical width of the target beam. 



3 Data Acquisition and Processing 

3.1 Measuring Cycle 

Measurements of the excitation functions A^ipiab, 0c.m) were 
performed in cycles of about 15 s duration with data acquisi- 
tion extending over the acceleration, the flattop at 3.3 GeV/c, 
and the deceleration as well. A typical machine cycle is shown 
in Fig. The beam current is nearly constant during the ac- 
celeration and slowly decreasing during the deceleration. The 
trigger rate and the dead time show huge excursions at certain 
COSY-times which are caused by increased background due 
to beam losses. Fortunately, these excursions occured at differ- 
ent COSY-times in the three running periods. The direction of 
the target polarization was changed cyclewise yielding the se- 
quence +x, —X, +y, —y. The target was operating very stable 
and with constant polarization during subsequent acceleration 
cycles. 

3.2 Identification of Elastic pp Events 

The outer detector shell provides a fast and efficient trigger 
based on (i) the coplanarity and (ii) the kinematic correlation 



1 




Fig. 8. Definition of the kinematic deficit a in the cm. system. 



of two-prong events that fulfill the kinematics of elastic proton- 
proton scattering. Elastic proton-proton scattering events are 
identified by coplanarity with the beam axis 



\ipi-^2\ - 180° 



and kinematic correlation 



tan 01 tan 6*2 = 



2mp + Tiab 



(2) 



(3) 



Here, and ip are the polar and azimuthal angles in the lab 
system, nip is the mass of the proton and Tiab its laboratory 
kinetic energy. In contrast to measurements of da/dil with CH2 
fiber targets ^ ^ the data rate with an atomic beam target 
was rather low. Therefore, only the coplanar trigger, i.e. the 
coincidence of two opposite bars, was applied in the on-line 
trigger for the analyzing power measurements. The kinematic 
correlation was established in the off-line analysis. 

The signature of an elastic event can be represented by one 
variable, the so-called kinematic deficit a, which gives the spa- 
tial angle deviation from back-to-back scattering in the c.m.- 
system (see Fig.|8}. The kinematic deficit a can be determined 
in the off-line analysis by transforming the trajectories of two 
coincident particles into the cm. -system assuming the kine- 
matics of elastic proton-proton scattering. The resulting dis- 
tributions of the spatial angle a start with zero at a = 0° and 
show a narrow peak followed by a long tail, see Fig.|9] The fi- 
nite width of the elastic peak is due to the effects of small angle 
scattering and the finite angle resolution of the detector. Elas- 
tic pp scattering events can be identified using a momentum 
dependent cut, 

a<aeut, ac„t = [11.0- p/(lGeV/c)]°. (4) 



3.3 Vertex Reconstruction 

In the off-line analysis the coordinates of two kinematically 
correlated trajectories are deduced from the hit pattern in the 
inner and outer detector shells. The vertex is determined ge- 
ometrically by the subroutine FINDTRACKS as the point of 
closest approach of two outgoing trajectories in the target re- 
gion. 

In order to improve the vertex reconstruction the data are 
reanalyzed in a kinematic vertex fit with the kinematic con- 
straints of elastic proton-proton scattering. The kinematic fit of 
two proton trajectories to the hit pattern of an event under the 
constraint of elastic scattering kinematics is used to define a 
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Fig. 9. Measured distributions of the kinematic deficit a at two beam 
momenta. The vertical lines show the momentum dependent cut for 
the selection of elastic pp scattering events. 



X criterion for a further event selection. This method yields 



for instance at r;af,=1500 MeV and 0c 



=40° a vertex reso- 
lution with ax=0-6 mm, ay— 0.6 mm and az—OA mm and an 
angle resolution with ctq^ ^ =0.6° and a^p— 0.8°. 

The x^-distribution is equivalent to the a-distribution with 
respect to the identification of elastic scattering events. The up- 
per part of Fig.[TO|shows an example of a x^-distribution from a 
kinematic vertex fit at p;a&=2280 MeV/c and 6'c.m.=89°. The 
dashed line is an extrapolation in order to estimate an upper 
limit of the background contribution. A momentum dependent 
cut 



2 , 2 



xLt - 28.0 - 5.5p/GeV/c 



(5) 



was chosen as final selection criterion for an elastic pp scatter- 
ing event, see vertical lines in Fig. 1101 



3.4 Background 

The reduction of background takes advantage of the recon- 
structed vertex and the multiplicity patterns in both detector 
layers. NaiTow cuts were applied to the hit pattern and to the 
vertex coordinate z in COSY beam direction, [-15 mm, +20 
mm], a wider one in the x, y plane around the beam profile 
(3(7 limits) along an ellipse following the slow drift and shape 
variation of the COSY beam during acceleration, see Figs. 15161 
The remaining inelastic background was estimated guided by 
Monte Carlo simulations of elastic and inelastic pp interac- 
tions. For the simulation of the hadronic reactions the code MI- 
CRES 1 67 1 was used. 

Monte Carlo simulations show that the long tail of the a- 
and -distributions is mainly due to misidentified elastic scat- 
tering events suffering from a secondary reaction in the beam 
pipe and the inner detector shell. This is in accordance with 
the fact that events with a>acut (x^>Xrat) show analyzing 
powers very similar to the elastic scattering. Therefore, extrap- 
olating the a- and x^-distributions from large x^-values to zero 
as shown in Fig.Qoloverestimates the inelastic background un- 
der the elastic peak. This background of inelastic reactions like 
pp — + ppT^^, PP —>■ pmr^, pp —^ ppn^TT^ , pp — > ppn'^Tr'^, and 
pp prni^Ti'^ is rather small. It was estimated to be mostly < 
2% and only at highest energies near 6>c.m.=90° up to 4.5%. 





Pl^b = 2280 MeV/c 
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Fig. 10. Upper part; Example of a x^-distribution. The dashed line 
is an extrapolation from x^>30 to zero. The vertical line shows the 
cut for the selection of elastic pp scattering events. Lower part: Monte 
Carlo simulation of a x'^ -distribution at 2.3 GeV/c. (a) Extrapolation 
from x'^>30 to zero, (b) The inelastic background contribution (lower 
histogram) 



3.5 Determination of Analyzing Power 

We denote the target polarization by Q={Qx, Qy,Qz)- In or- 
der to eliminate systematic errors the direction of the target po- 
larization was changed from cycle to cycle between ±x and 



Q\ 
Q\ 

y - 
Qy=-\Q\- 

The polarized differential cross section may be written 



+ x 


Q = 


(+ 


QI,o,o), 


— X 




(- 


QI,o,o), 


+y 


Q = 


(0, 


+IQI,o), 


-y 


Q = 


(0, 


-|QI,o), 



A-N {Qy COS ip- Qx sin tp)] . 



(6) 



(7) 



Here, dco/d]? is the unpolarized differential cross section. An 
the analyzing power with respect to a polarization component 
in the direction N normal to the scattering plane. The azimuthal 
angle tp gives the rotation of the scattering plane and it's coor- 
dinate system (5, N, L) with respect to the fixed coordinate 
system (x,y,z) of beam and target, see Fig. II II 

In order to eliminate false asymmetries arising from differ- 
ences of the luminosities, efficiencies and solid angles of the 
detector and from misalignments, the geometric mean method 
of Ohlsen and Keaton ll63ll is used. Two sets of cycles with op- 
posite polarizations, e.g. (5+ and Q", were combined to apply 
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Fig. 11. Definition of four detector sectors for the determination 
of the analyzing power with, looking along the beam axis, L (Left), 
R (Right), B (Bottom) and T (Top). The azimuthal angle if gives the 
rotation of the scattering plane and it's coordinate system {S, N, L) 
with respect to the fixed coordinate system (i, y, z) of beam and tar- 
get. The sectors with 75°<i/j<105° and 255°<ip<285° are not used 
since the angle reconstruction is affected by the readout of the half 
rings, see Fig. 1. 



a "proper spin flip" correction for false asymmetries. Indicating 
the number of events obtained simultaneously for the detector 
elements in the angle position (6*, ip) with N^{0, (p) for spin 
up and {0, (p) for spin down we define 



L = y^N+{0,ip)N-{0,ip + Tr) 



and 



i?= ^N-{e,ip)N+{0,ip + 7r) 
to calculate the left-right asymmetry 

L- R 



L + R 



(8) 



(9) 



(10) 



for a pair of detector elements in the azimuthal positions ip and 
V? + TT with — 7r/2 < ip < +tt/2. The analyzing power An{0) 
is deduced as weighted mean over all p-hins using 



ANiO) = 



1 ^LB.{6>,ip) 



COS ip 



Qy 



^7r/2 < (/? < 7r/2. (11) 



Here, the moduli of the target polarizations are assumed to be 
equal,|Q+| = |Q-| = |Q|. 



Similarly the runs with were used to deduce An from 
the bottom-top asymmetries esT for a pair of detector ele- 
ments in the azimuthal positions (p and (p+ir with — tt < (p < Q, 



ANiO) 



I eBT{0,v) 
sin Lp Qx 



-TT < < 0, 



assuming |Q+| = |g-| = |Q|. The terms Left (L), Right (R), 
Bottom (B) and Top (T) always refer to the scattered proton 
detected at forward scattering angles with Oc.m. < 90°. 
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Fig. 12. Angular distribution at 1379 MeV/c (730 MeV) in compari- 
son to the LAMPF data J5§\ which were used for the absolute normal- 
ization. The momentum bin width of the EDDA data is here 60 MeV/c. 



This method eliminates exactly all false asymmetries that 
means asymmetries which would still be observed with no tar- 
get polarization. Thus, the result is independent of relative de- 
tector efficiencies and solid angles, since they do not vary with 
time over the period of two adjacent cycles. It is also indepen- 
dent of time fluctuations in the beam current or target density 
as well as differences in the integrated charge and target thick- 
ness. 

Misalignments of the beam axis with respect to the detector 
axis yield small deviations from the nominal scattering angles 
Oc.m. and small variations of the solid angles. These misalign- 
ments depend on the closed orbit distortions during the beam 
acceleration and deceleration. Again, there is exact cancellation 
of false asymmetries due to deviations of the solid angles. But 
since the analyzing power An depends on the scattering angle 
Oc.m. the determination of the effective scattering angle be- 
comes an important experimental task. Fortunately, the EDDA 
detector allows to reconstruct the scattering angle for each 
event with high accuracy. Thus, systematic errors from these 
deviations are also avoided. 

The geometric mean correction method assumes the mod- 
uli of the target polarizations and to be equal. This 
assumption is very well fulfilled since the polarization vector 
follows adiabatically the direction of the magnetic guide field 
and the spin flip is realized by flipping the direction of the guide 
field. Small deviations that may occur cause negligible effects. 
For instance deviations with \ — \Qy \ < 0.02 influence An 
by at most 7 • 10~^. 



3.6 Absolute Target Polarization 

The absolute values of the effective target polarizations and 
Qy are established in each running period for one momentum 
bin Apiab = 60 MeV/c around the energy Tiab = 730 MeV 
(12) (Piab — 1379 MeV/c). The absolute values of Qx and Qy are 
deduced as weighted means over all (/j-bins from the measured 
asymmetries. 



Q. 



sin Lp An (O) ' 



-TT < ip <0, 



(13) 
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Qy 



-n/2 < ip < 7r/2. 



(14) 



cos (pAn{0) ' 

The precise angular distribution An{0) from McNaughton et 
al. I55I1 with an absolute normalization error of 1 % is taken as 
reference value, see Fig. [21 Taking the relative errors of the 
data into account the overall normalization error of the present 
data is 1.2 %. The effective target polarization is constant dur- 
ing the acceleration and deceleration since the overlap between 
COSY beam and polarized atomic beam target is constant. The 
observed small variations of the vertical beam width (see Fig.|6j 
cause variations of the effective target polarization of less than 
0.3 %. 




tODO t250 tSOO tT50 2000 2250 2500 2750 3000 3250 

Pi*(MeV/c) 



3.7 Effective Target Polarization at Forward Angles 

The influence of the restrictions by the detector acceptance has 
been studied separately. It yields small modifications of the ef- 
fective target polarization at low energies and small scattering 
angles. There, events with negative z-values of the vertex point 
cannot be detected if the recoil protons are outside the detec- 
tor acceptance. As a consequence the effective target volume 
is reduced in the z-direction from [—15 mm, +20 mm] to 
[zmin, +20 mm] and the average polarization within this re- 
duced volume, subsequently called the effective polarization, 
differs from the full-acceptance value. This modification is due 
to the background of unpolarized hydrogen atoms, see Fig.|3] It 
can be studied in a systematic way by artificially restricting the 
vertex range in z-direction for all {Oc.m., Piab)-^ins with full 
acceptance and deducing the weighted mean of the resulting 
asymmetry ratios e{[zrain, +20 mm])/e([— 15 mm, +20 mm]) 
as a function of z,nin- These asymmetry ratios can directly be 
used as correction factors of the effective target polarization at 
low energies and small scattering angles where the z-range is 
restricted. 



4 Results 

4.1 Consistency Checks 

The COSY beam changes its shape and position during accel- 
eration and deceleration, though very reproducible in each cy- 
cle. Prior to merging all data in one final set it was necessary 
to perform consistency checks on subsets obtained under dif- 
ferent conditions of polarization and acceleration cycle. They 
demonstrated that the guide field is properly aligned to the de- 
tector coordinates (x,y), and that during acceleration and de- 
celeration the same analyzing powers are obtained |2]. This 
implies that vertex reconstruction and proper flip elimination 
of false asymmetries work well. A comparison of analyzing 
powers An at 0c.m.—GQ° acquired during the running periods 
June 98, November 98 and May 99 is shown in Fig. ^] Since 
all data are compatible they are combined in one set. Altogether 
3.1-10^ elastic scattering events were collected. 

4.2 Errors 

Error estimates for An include contributions from the maxi- 
mum deviations between data subsets (< 2.1%), and the im- 
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Fig. 13. Comparison of analyzing powers ^jv at 0c.m.~6O° acquired 
during the running periods June 98, November 98 and May 99. 



pact of the background on the asymmetry (< 0.008). Errors 
due to closed orbit distortions by the variation of the magnetic 
guide field with changes of the proton beam position and angle 
are negligible, see Sect. 2.4. The overall absolute normaliza- 
tion uncertainty of the excitation functions is ±1.2 % in the 
full energy range 0.45-2.5 GeV. It is due to the absolute nor- 
malization uncertainty of ±1 % from the reference data Issll 
and the relative errors of the data at 730 MeV. This systematic 
uncertainty is not included in Figs. 1141 and[T6land the tables 
available upon request itssil . and must be apphed to all data. 



4.3 Excitation Functions 



Excitation functions A^ipiab, ©cm.) with about 1 150 data have 
been deduced from our experimental results by grouping them 
into A0c.m.^^° and Apiab—'iO MeV/c wide bins. They super- 
sede the results of the EDDA Collaboration reported in |2i| and 
are available upon request l'68'l. Here, excitation functions at 
eight out of 15 cm. angles are displayed in Fig. 1141 



4.4 Angular Distributions 



The data can also be presented in the form of angular distri- 
butions. In Fig. eight out of 77 angular distributions for 
A0c.m.~^° and Apiab—iQ MeV/c wide bins are shown. Again, 
these data supersede the results of the EDDA Collaboration re- 
ported in |2|] and are available upon request 16 8r . 
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Fig. 14. Excitation functions of An from the present work at eight Fig. 15. Collection of published data from the data base of fl^ with- 
out of 15 cm. angles in comparison to the phase shift solutions out EDDA data of Altmeier et al. |2] plotted as excitation func- 
SPOO (dashed curve) and FAOO (solid curve) |Tj]. The bin widths are tions at eight cm. angles in comparison to the phase shift solutions 
Apiab=SO MeV/c and A0c.m.='i° ■ The data at the flattop momen- SPOO (dashed curve) and FAOO (solid curve) fl3|. O SATURNE II 
tum 3.3 GeV/c exhibit an especially small statistical error. IHIilE^lsSEIlIm ^ KEK 1 61], O ZGS 1 69jc3l2ill. ^ LAMPF 

mUlIll, □ lUCFBai, + SIN |73|73H7J. 



5 Discussion 

5.1 Comparison to other data 

Previous An data were mainly measured at discrete energies. 
A collection of published data from the SAID data base is 
shown in Figs.^]and[n] The most recent results in the kinetic 
enemy range 0.8-2.8 GeV are from SATURNE II S El US, 
l5lll52il . They are in good agreement with the present data. The 



same holds true for the results from SATURNE II of Perrot et 
al. |47] at the kinetic energies (momenta) 0.874 (1.550), 1.095 
( 1 .804), 1 .295 (2.027), 1 .596 (2.354), 1 .796 (2.568), 2.096 (2.886) 
and 2.396 GeV (3.200 GeV/c). The ZGS data of Parry et al. 
ii at 1 .73 (2.5), 1 .97 (2.75), 2. 14 (2.93) and 2.44 GeV (3.25 GeV/c) 
and Miller et al. jH at 1.27 (2.0) and 2.21 GeV (3.0 GeV/c) 
show considerable deviations in the angular distributions. 
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Fig. 16. Angular distributions of Am from the present work at eight 
out of 77 beam momenta in comparison to t he p hase shift solutions 
SPOO (dashed curve) and FAOO (solid curve) Il4r . The bin widths are 
Apiab=30 MeV/c and A0c.m.=^° ■ The data at the flattop momen- 
tum 3.3 GeV/c exhibit a small statistical error. 



Fig. 17. Collection of published data from the data base of with- 
out EDDA data of Altmeier et al. |2] plotted as angular distribu- 
tions at eight momenta in comparison to the phase shift solutions 
SPOO_(dashed curve) and FAOO (solid curve) flil O SATURNE II 
IHrarais^EllsiLA KEK|61],OZGS 1 69^03 IZill-^LAMFF 
milllll, □ lUCFlHi, + SIN |73|73H7J. 



In the kinetic energy range 0.45-0.8 GeV the world data 
set exhibits several high precision measurements at discrete en- 
ergies. The absolute normalization of our An excitation func- 
tions was established at r/Qb=730 MeV (p/Qb=1379MeV/c) by 
taking the precise angular distribution An{0) of McNaughton 
et al. [5.51 as reference value. Our data below and above 730 
MeV (1379 MeV/c) are in very good agreement with other 
precise LAMPF measurements of McNaughton et al. ll54ll at 



643 MeV (1273 MeV/c) and Bevington et al. (sf] at 787 (1448) 
and 796 MeV (1459 MeV/c) as well as with precise lUCF mea- 
surements of Przewoski et al. |58| at 448.9 MeV (1022 MeV/c), 
SIN measurements of Besset et al. |72], Berdoz et al. 1 .7 3ll . 
Aprile et al. |74| between 400 and 600 MeV (954 and 1219 
MeV/c) and SATURNE measurements of Allgower et al. (mI 
at 795 MeV (1457 MeV/c). 
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5.2 Comparison to Partial-Wave Analyses 

After the publication of the first unpolarized differential cross 
section data from EDDA (l| the VPI group extended their en- 
ergy dependent phase shift analysis from 1 .6 GeV (2.36 GeV/c) 
up to 2.5 GeV (3.3 GeV/c) laboratory kinetic energy (beam 
momentum) with the solution SM97 |13]. Meanwhile energy 
dependent phase shift solutions are available with a maximum 
beam energy (beam momentum) of 3.0 GeV (3.82 GeV/c) 1 14]. 
The solution SPOO includes the recent data from lUCF 
156, .57. 58., .59.1 and Saclay |^^. The recent An data from 
EDDA l^are included in the solution FAOO. 

The comparison of the present data to the phase shift so- 
lution FAOO 1 14] yields agreement in the size and the general 
angle and momentum dependence of the excitation functions 
and angular distributions, see Figs. 1141 andll6l Small system- 
atic deviations can be seen in the excitation functions, in par- 
ticular for momenta from 1.8 - 2.5 GeV/c. As can be seen in 
Fig. the difference between the solutions FAOO and SPOO 
is small and the systematic deviations in the momentum range 
1800 - 2500 MeV/c remain. It is interesting to note that the 
analyzing power data are slightly negative in the region around 
Pzab=2.0 GeV/c and ©cm. =80° whereas the the phase shift so- 
lutions remain positive. This observation is in agreement with 
previous Saclay data |47] in that region. 

Including the EDDA data into the solution FAOO turned out 
that the agreement of the phase shift solution with the new A n 
data was slightly improved. The spin triplet phases, e.g. those 
of the partial wave experienced significant changes. This is 
due to the fact, that the analyzing power times differential cross 
section is equal to the real part Re{a* ■ e), where the invariant 
amplitudes a and e 0,0] include only triplet partial waves. 



5.3 Sensitivity to narrow resonances 

All excitation functions of the analyzing powers show a smooth 
dependence on beam momentum. There is no evidence for nar- 
row resonances. 

A previous internal target experiment at KEK 1 6 1 ] observed 
two narrow structures in the excitation function of the analyz- 
ing power An at a laboratory kinetic energy near 632 MeV 
corresponding to -^s = 2.17 GeV. Those measurements were 
performed using a polarized proton beam in the KEK ring and 
a 30 /Ltm thick polyethylene fiber target. The outgoing pro- 
tons were detected at a fixed backward angle of 68° in coin- 
cidence with a forward detector The corresponding cm. angle 
near 632 MeV was 38.4°. In the present work all excitation 
functions including those near 38° are smooth. Our result is in 
agreement with the measurements of Beurty et al. i62ll . 

The fact that ANda/dfl = Re{a* ■ e) includes only triplet 
partial waves implies that the excitation functions for An are 
more sensitive to resonant excursions in triplet than to those in 
singlet partial waves. However, the excitation functions mea- 
sured here exhibit no evidence for energy-dependent structures. 
It should be noted that also the excitation functions of da/dfl 
did not show any evidence for energy-dependent structures 
Qi.A more detailed discussion of sensitivities to and upper 
limits for such structures will be given in a forthcoming paper. 



6 Conclusions 

In conclusion, we report on the first measurement of analyz- 
ing power excitation functions AN{piab, Oc.m.) in the labora- 
tory momentum range 1.0-3.3 GeV/c and the cm. angle range 
30°-90° for proton-proton scattering during acceleration and 
deceleration in a synchrotron. The data provide a new polariza- 
tion standard and can be used for calibration purposes in the 
full energy range 0.45-2.5 GeV. The excitation functions agree 
with fixed energy data and close the gaps in between with data 
of high precision and consistency. The phase shift analysis in- 
cluding our data yields a global phase shift solution for Tiab 
up to 2.5 GeV (FAOO) showing distinct deviations from previ- 
ous phase shift solutions that occur mainly in the spin triplet 
phases. Further progress can be expected from the new exci- 
tation functions of the spin correlation parameters Ann, Ass 
and AsL measured by the EDDA experiment. 
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